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Abstract 

In this paper, we consider a system of k second order non-linear stochastic differential 
equations with spatial dimension d > 1, driven by a fc-dimensional Gaussian noise, 
which is white in time and with some spatially homogeneous covariance. We prove 
existence, smoothness, and strict positivity of the density of the law of the solution 
of this system of equations, on the set where the diffusion matrix is invertible, under 
sufficient conditions on the fundamental solution T of the deterministic equation. For 
this, we apply techniques of Malliavin calculus. We apply this result to the case of the 
stochastic heat equation in any space dimension and the the stochastic wave equation 
in dimension d 6 {1,2,3}, with a spatial covariance given by a Riesz kernel. We then 
study the strict positivity of the density for the case of a single equation (k = 1), and 
apply it to the stochastic heat equation in any space dimension, and the stochastic wave 
equation in dimension d s {1, 2, 3}, with a general spatial covariance. 
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1 Introduction and main results 



Consider the following system of stochastic partial differential equations: 

k 

Lut(t,x) =^2a ij (u(t,x))W k (t,x) + bi(u(t,x)), t > 0, x G R d , (1.1) 

3=1 

i = 1, k, with initial conditions 

Ui (0,x) = ^(0,x) = 0. 

Here L is a second order differential operator and <7y , : M fe i— ► R are globally Lipschitz 
functions, that is, there exists a constant K > such that for all x,y £ W 1 and 1 < < fc, 

| cry (a;) - o-y (y)| + - < K\\x - y\\. 

We denote by a := ((Tij)i<i,j<k the diffusion matrix and 6 = (6j)i<j<fc the drift vector. The 
driving perturbation W"(i, x) = (W^(t,x), ...,W k (t,x)) is a /c-dimensional Gaussian noise 
which is white in time and with a spatially homogeneous covariance /, that is, 

E [W%x)Wi(s,y)] = 5(t - s)f(x - y)S ij , 

where <5(-) denotes the Dirac delta function, 5ij the Kronecker symbol, and / is a non- 
negative continuous function of M. d \ {0} such that it is the Fourier transform of a non- 
negative tempered measure jj. on R rf . That is, 



f(x) = ^^x):= [ e-^ M (de). 



Moreover, we assume that for some integer m > 1, 



12^™ 



< +00. 



Let {^t)t>o denote the filtration generated by W (see Section 2 for the precise definition 
of the noise), and let T > be fixed. By definition, the solution to the formal equation 
(|l-ip is an adapted stochastic process {u{t,x) = (ui{t,x),...,u k (t,x)),(t,x) G [0,T] xR d } 
such that 

Ui(t,x) = y2 / r (* - s ' x - y)o-ij{u{s,y))W 3 {ds,dy) 

3=1 JoJ ** (1.2) 

+ / / bi(u(t - s,x - y))T(s,dy)ds, 
Jo JR d 

where T denotes the fundamental solution of the deterministic equation Lu = 0. 

Recall that when T(t,x) is a real valued function, the stochastic integral appearing in 
(|1.2|) is the classical Walsh stochastic integral (see W86J). However, when T is measure, 
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Dalang [D99] extended Walsh's stochastic integral using techniques of Fourier analysis, and 
covered, for instance, the case of the wave equation in dimension three. D.Nualart and 
Quer-Sardanyons [NQ07] extend Walsh's stochastic integral using techniques of stochastic 
integration with respect to a cylindrical Wiener process (see [DZ92] ). in order to cover 
the case of measure- valued integrands. This theory will be recalled in Section 2, in our 
/c-dimensional setting. This extension turns out to be equivalent to Dalang's integral when 
the integrands are of the form G := T(t — , x — *)Z(-, *), for certain stochastic processes Z 
(see Section 2). 

We now introduce some hypothesis on T, and on the coefficients of equation (|1.2p . that 
wil be needed along this paper. 

(HI) For all t > 0, T is a nonnegative distribution with rapid decrease, such that 

f [ \&T(t)($)\*ii(d£)dt <+oo. (1.3) 
Moreover, T is a nonnegative measure of the form T(t, dx)dt such that 

sup / T(t,dx) < C T < +oo. (1.4) 

te[o,T] JR d 

(H2) The functions Uij and bi are < ^' 00 bounded functions with bounded partial derivatives. 
(H3) There exists 71 > such that for all h £ [0, T], 

[ [ \&Y{h + r){£) -&T(r)(£)\ 2 ii(d£)d,r < c T h 71 , (1.5) 

Jo JR d 

for some positive constant ct > 0. There exists 72 > such that for all h G [0, T], 

f [ \&T(r)(i)\ 2 ix{d£)dr <c T h?\ (1.6) 
Jo JR d 

for some positive constant ct > 0. There exists 73 > such that for all h € [0,T — t], 

r-t+h f 

(1.7) 



rt+n p 

/ / T{s,dy)ds < c T h ri , 

Jt JR d 



for some positive constant ct > 0. Set 7 := 71 A 72 A 273. 
(H4) There exists 74 > such that for all y, z G M. d , 

f T [ \#T(r)(y -0- &T(r)(z - i)\ 2 ^d£)dr < c T \\y - z\ 
JO JR d 

for some positive constant ct > 0. 



74 



3 



(H5) The nonnegative measure ^ defined as ||x|| 2 T(t,dx)dt satisfies that 



/ / \&9(t)(0\* < +00. 



JO JR d 



Moreover, there exist a.\,ot2 > 0, ^ < a\, % < 02, with ot\ A ai < (272) A (72 + 1), 
such that for all r € [0, 1], 



for some positive constant c > (see Section 2 for the definition of the space and 



for some positive constant c > 0. 

In the paper by D.Nualart and Quer-Sardanyons N Q"07] , the authors use techniques of 
Malliavin calculus in order to obtain existence and smoothness of the density of the law 
of the solution to (jl.ip in the case where k = 1. For this, they assume hypotheses (HI), 
(H2), (H6), and the fact that |er| > c> 0. 

The first aim of this paper is to prove existence and smoothness of the density of the law 
of the solution to (II. ip under conditions (Hi), i = 1, 6, on the set where a is invertible, 
using techniques of Malliavin calculus (see Theorem 11.1( a) ) . This will be done in Section 
4. For this, we use a natural fc-dimensional extension of the results in |NQ07|, together 
with a localizing argument similar as the one used in |BP 98] . This requires the additional 
assumptions (H3), (H4), and (H5). We observe that if we assume that the matrix a is 
uniformly elliptic, one also obtains the uniformly boundedness of the density (see Remark 



Observe that hypotheses (H3) and (H4) imply in particular that the paths of the 
solution to (jl.ip are 7/2-Holder continuous in time, and 74/2-Holder continuous in space 
(see Lemma I3.2p . The Holder continuity of the paths of solutions of equations of the type 
(jl.ip have been largely studied in the literature. The case of the stochastic heat in any 
space dimension was studied in [MS99| . and the stochastic wave equation in two space 
dimension was developed in |SS02j . For Hdlder-Sobolev regularity results for the stochastic 
wave equation in three space dimension see [DS09J. 

The second and principal aim of this article is to show that under the same conditions 
(Hi), i = 1,...,6, the density of the law of equation (jl.ip is strictly positive on the set 
where a is invertible (see Theorem 11.1( b)). For this, we will first extend in our situation a 
criterium of strict positiveness of densities proved by Bally and Pardoux in |BP98j . which 





for some positive constant c > 0. Set a := a\ A 02- 



(H6) There exists r\ > with < 72 < rj < a such that for all r G [0, 1] 




El. 
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uses essentially an inverse function type result (see Lemma I5.1|) and Girsanov's theorem. 
This will be done in Section 5. 

Section 6 will be devoted to apply the criterium of strict positiveness of Section 5 to our 
class of systems of SPDEs (jl.ip under the mentioned conditions. In [BP98J, the authors 
apply their criterium to the density of the law of the random vector (u(t, xi), ....,u(t,X)-)), 
< X\ < • • • < Xk < 1, where u denotes the solution to the non-linear stochastic heat 
equation driven by a space-time white noise, by using a localizing argument. Hence, their 
situation is a bit different as ours, as we deal with a system of SPDEs and we evaluate the 
solution at a single point (t, x) G]0, T] x M. d . In a similar context, Chaleyat and Sanz-Sole in 
[CS03J study the strict positivity of the density of the random vector (u(i, xi), u(t, x^.)), 
< x\ < • • • < Xk < 1, where u is the solution to the stochastic wave equation in two 
spatial dimensions, that is, take in equation (jl.ip d = 2, k = 1, and L the fundamental 
solution of the deterministic wave equation. 

We will also apply the criterium of strict positivity of the density of Section 5 to the case 
of a single equation, that is the solution to (jl.ip with k = 1. As we have recalled above, in 
[NQ07], the authors prove existence and smoothness of the density of this equation under 
hypotheses (HI), (H2), (H6), and the fact that |<t| > c > 0. We prove that the density is 
strictly positive in all M. under the additional hypothesis (II. 6p (see Theorem 11.31 proved in 
Section 7). 

In Section 8, we will apply our results to a system of stochastic heat and wave equations 
driven by a spatial covariance defined as a Riesz kernel (see Theorems 18.11 and I8.3P . Observe 
that in the wave equation case, only in the case tZ G {1,2,3} T defines a nonnegative measure. 
We also study the case of a single stochastic heat equation in any space dimension (see 
Theorem 18. 2p . and a single stochastic wave equation in spatial dimension d € {1,2,3} (see 
Theorem I8.4D . all driven by a general spatial covariance. 

There are also other many SPDEs for which the strict positivity of its density is studied. 
For example, the case of non-linear hyperbolic SPDEs has been studied by Millet and Sanz- 
Sole in |MS97j . Fournier [F99j considers a Poisson driven SPDE, and the Cahn-Hilliard 
stochastic equation is studied by Cardon- Weber in [C02J. 

One of the motivations of the results of this article, is to develop in a further work 
potential theory for solutions of systems of the type (II. ip (see [DKN09b|). For the moment, 
potential theory for systems of non-linear SPDEs has been studied by Dalang and E.Nualart 
in [DN04] for the wave equation, and by Dalang, Khoshnevisan and E.Nualart in [DKN08J 
and [DKN09a], for the heat equation, all driven by space-time white noise. That is, taking 
in equation (jl.ip . d = 1, / the Dirac delta function and L the fundamental solution of 
the deterministic wave or heat equation. In all these works, the existence, smoothness, 
uniformly boundedness, and strict positivity of the density of the solution to the system of 
SPDEs is required. 

We next state the main result of this paper: 

Theorem 1.1. Assume hypotheses (Hi), i = 1,...,6. Then for all (t,x) €]0,T] x W 1 , the 
law of random vector u(t,x) admits a c /a QC strict positive density on the open subset ofM. k 
£ := {y G R k : det a(y) / 0}. That is, 

(a) there exists a function pt, x € ^^(SjR) such that for every bounded and continuous 
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function f : M. k i— > R ot^/j support contained in E, 

E[/(«(i,x))] = / f(y)pt^{v)dy, 

(b) and pt,ar(j/) > for all y G E. 

Remark 1.2. If a is strongly elliptic, that is, • £|| 2 > p 2 > (or, equivalently, since 

a is a square matrix, ||£ T • a(x) \\ 2 > p 2 > 0) for some p > 0, for all x G R d , £ G R d , ||£|| = 1, 
then one obtains the additional property that the density pt,x is uniformly bounded in any 
non-degenerate compact rectangle Ax B c]0, T] x W 1 . 
Note that because a is a square matrix, it holds that 

inf inf \\a(x) ■ £|| 2 = inf inf ||£ T • a(x)\\ 2 . 
xeK d ||5||=i xeR d ||5||=i 

However, for non square matrices this equality is false in general. 

We finally state the case of a single equation {k = 1). 

Theorem 1.3. Let {u(t,x),(t,x) € [0, T] x M. d } denote the solution of equation with 
k = 1. Assume hypotheses (HI), (H2), (II, 6p . and (H6). Suppose also that \o~(z)\ > c > 0, 
/or a// z£l, Taen /or a// (t, x) G]0,T] x the law of random variable u(t,x) admits a 
strictly positive density on R. 

2 The stochastic integral 

The aim of this section is to recall the results obtained in N Q07| in our fe-dimensional 
situation. Fix a time interval [0, T]. The Gaussian random perturbation W = (W 1 , W k ) 
is a fc-dimensional zero mean Gaussian family of random variables W = {W J ((/?), 1 < j < 

^G^dO.Tjxl^R)}, 

The aim of this section is to recall the results obtained in |NQ07| in our fc-dimensional 
situation. Fix a time interval [0, T]. The Gaussian random perturbation W = (W 1 , W k ) 
is a fc-dimensional zero mean Gaussian family of random variables W = {W^p), 1 < j < 
d, ip G < ^ ,o ([0,T] x B^jM)}, defined on a complete probability (O, J^,P) with covariance 

E[r(^)^W] = 5ij f T j j <p(t,x)f(x - z)il>(t,y) dxdydt. 

JO JR d JR d 



This covariance can also be written as 



T 







Let Jt? k denote the completion of the Schwartz space y(R d ; R k ) of <?f 00 (R d ; R fc ) functions 
with rapid decrease, endowed with the inner product 

M)** = J2 [ I Mx)f{x - y)Mv) dxdy = V / &M0&MtM<%), 

Tl J^ d JR d e, JR d 
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4>,i> € y(R d ;R k ). Notice that may contain distributions. Set Jt? T = L 2 ([0, T];J4? k ). 
In particular, Ji?r is the completion of the space ^^([O, T] x R d ;R fc ) with respect to the 
scalar product 

M)^, = V [ T [ & n {tW)^Mm)m)dt. 

Then the Gaussian family W can be extended to Mr and we will use the same notation 
{W(g),g G J? T }, where W(g) = £* = i W l ( gi ). 

We next recall the extension of Walsh's integral defined in |NQ07| . Set Wt(h) = 
EiLiW^lp^/i;), for any t > 0, h G J^f k . Then {W t ,t G [0,T]} is a cylindrical Wiener 
process in the Hilbert space J4? k (cf. [DZ92} Section 4.3.1]). That is, for any h G J4? k , 
{Wt(h),t G [0, T]} is a Brownian motion with variance \\h\\^ k , and 

E[W t (h)W t (g)] = (sAt)(<p,i>)^ k . 

Let (^t)t>o denote the cr-field generated by the random variables {W s (h), h G J4? k , < s < 
i} and the P-null sets. We define the predictable cr-field as the cr-field in Q x [0, T] generated 
by the sets {(s,t] x A,0 < s < t < T, A £ & s }. Then following |DZ92l Chapter 4], we can 
define the stochastic integral of any predictable process g G L 2 (Q x [0, T];Jif k ) with respect 
to the cylindrical Wiener process W and we denote it by 



g ■ dW = g ■ W. 
Moreover, the following isometry property holds: 









' r 2 


E 


\g-W\ 2 


= E 


/ \\9t\\^kdt 








Jo 



Then using |NQ07 Proposition 3.3], one can define the stochastic integral of a predictable 
process G = G(t,dx) = Z(t, x)T(t, dx) G L 2 (Q x [0,T]; Jf k ) with respect to W, provided 
that Z = {Z(t,x), (t,x) G [0,T] x R. d } is a predictable process such that 

sup E[\\Z(t,x)\\ 2 ] < +oo, 

(t,x)e[0,T]xR d 

and in this case, we denoted the stochastic integral as 



G-W 



G(s,y)-W(ds,dy) 



T(s,y)Z(s,y)-W(ds,dy). 



3 Existence, uniqueness and Holder regularity of the solution 

We assume that for all i = l,...,k fixed, the process Zj(s,y) = aij(u(s,y)), j = l,...,k 
satisfies the hypothesis of |NQ07 Proposition 3.3] so that the stochastic integral in (II. 2D is 
well-defined. 

We next state the existence and uniqueness result of the solution to equation (|1.2|) : 
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Theorem 3.1. NQ07, Theorem 4.1] Under condition (HI), there exists a unique adapted 
process {u(t, x) = (u\(t, x), Uf.{t, x)), (t, x) G [0,r]xR rf } solution of equation (|1,2|) . which 
is continuous in I? and satisfies that for all T > and p > 1, 

sup E[\\u(t,x)\\ p ] < +oo. (3.1) 

(t,x)£[0,T]xR d 

We next prove a consequence of hypotheses (H3) and (H4), which concerns an upper 
bound for the pth-moment of the increment of the process u, and in particular, gives the 
Holder continuity exponent of the paths of u. 

Lemma 3.2. Assume that condition (HI) is satisfied. 

(i) Under condition (H3), for allt E [0, T], h G [0, T — t], x € M. d , and p > 1, there exists 
7 > such that 

E[\\u{t + h,x)-u(t,x)\\ p ] <c p ^ K h lp/2 , (3.2) 

for some constant c P) t,k > 0. In particular, the trajectories of u are j/2-Hdlder 
continuous in time. 

(ii) Under condition (H4), for all t € [0, T], x,y G M. d , and p > 1, there exists 7' > 
such that 

E[\\u(t,x)-u(t,y)\\ p ] <c p ^ K \\x-y\\^l\ (3.3) 

for some constant c p ,t,k > 0. In particular, the trajectories of u are ^^jl-Holder 
continuous in space. 

Proof. We start proving (|3.2p . We write 

u(t + h, x) - u(t, x) = Ii + I 2 + h + h, 

where 

k ft p 

^ = J2 / (r(t + /i-s,x-y)-r(t- S ,x-y))(j,K S ,y))^(d S ,d2/), 

& pt+h p 

I * = J2 / r(t-s,x-y)a>(s,y))^(ds,dy), 

Jt JR d 

/ (6(u(t + h — s, x — y)) — b(u(t — s, x — y)))T(s, dy)ds, 
JR d 

t+h 

b(u(t — s, x — y))T(s, dy)ds. 

Using |NQ07 : , (3.11)], the Lipschitz property of the coefficients of <r, (|3.ip . and hypothesis 
(jl.5|) . we get that 



E[\\I x f]<c p , K [l+ sup E[\\u{s,y)f] 

(s,y)e[o,T]xR d 



T r \ p/2 

\&T{h + s)(0 - &T{s){t)\ 2 n{dt)ds 







< C v K T h^l\ 
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In the same way, but appealing to hypothesis (|1.6p . we obtain that 

E[||/ 2 H<c p Jl+ sup E[\\u(s,yW])( [ k [ \&T(s)(0\ 2 K^dsY^ 

V (s,y)e[0,T]xR d / \Jo JR d J 

< C p , K ,Th™ /2 . 

On the other hand, using the Lipschitz property of the coefficients of b, together with NQ07 
(5.17)], and hypothesis (jl.4|) . we get that 

EOlJslH <Cp, T ,K ! supE[\\u{t + h-s 1 y)-u(t-s,y)\\ p ] [ T(s,dy)ds 

JO y£R d JR d 

<c P)T ,K / sup E[\\u(h + s,y) - u(s,y)\\ p ]ds. 
Jo v m d 



Now, using the Lipschitz property of the coefficients of 6, together with Minkowski's inequal- 
ity with respect to the finite measure T(s, dy)ds, and appealing to (13. 1|) and hypothesis (|1.7p . 
we obtain that 

E[||/ 4 |n <c p Jl+ sup E[\\ U (s,y)\\ p ])( f +k [ T(s,dy)ds^' 

V (s,v)£\0,T]xR d / \Jt JR d 



< c p , K , T h™. 
Hence, we have proved that 



y)-u(s,y)\\ p ]ds). 



sup E [||«(* + h,y)- u(t, y) \\ p ] < c p , T ,i< (h? p ' 2 + f sup E[\\u(h + s, 

y£R d \ Jo y£R d 

Finally, apply Gronwall's Lemma to conclude the proof of (|3.2|) . 

We next prove (j3.3j) . Using |NQ07, (3.11)] and the Lipschitz property of the coefficients 
of a and 6, it yields that 



E[\\u(t,x)-u(t,y)\\P]<c p , K [l+ sup E[\\u(s,zW} 



(s,z)6[0,T]x] 



f T [ \^T(s)(x -0- ^E(s)(y - Z)\ 2 KdC)ds 

Jo JR d 



+ c p kE 



\u(s, x — z) — u(s, y — z) ||r(s, dz)ds 



p/2 



o JR d 



We now apply hypotheses (fPj) and (H4), and (pTTj) . together with |NQ07[ (5.17)], to get 
that 

E[\\u(t,x)-u(t,y)\\ p ] 

<c p , t ,k\\ x ~ y\V iP/2 + c p,t,k / sup E[\\u(s,x - z) - u(s,y - z)\\ p ] / T(s,dz)ds 

Jo z£R d JR d 

<c P)T ,k\\ x - y\\ lip/2 + c p,t,k / sup E[\\u(s,x - z) - u(s,y - z)\\ p ]ds. 

Jo z£R d 



In particular, this implies that 

sup E[||u(i 5 x — z) — u(t,y — z)\\ p ] 

zeR d 

rt 



+ c p ,t,k / sup E[\\u(s,x - z) - u(s, y - z) \\ p ]ds. 

Jo zm d 



< c p , t ,k\\x - y\r p / 2 
Finally, apply Gronwall's Lemma to conclude the proof of (13. 3p . 



□ 

We end this section proving a result which is a consequence of Lemma l3.2l and hypothesis 
(H5), that will be needed for the proof of Theorem 11.11 Its proof follows similar ideas as 
the ones in |NQ07, Proposition 3.3]. 

Lemma 3.3. Assume that hypotheses (HI), (H3), (H4), and (H5) are satisfied. Then 
for all t G [0, T], x € M. d , e G (0, 1], and p>l, 



E 



(aij(u(t - r, *)) - aij(u(t,x)))T(r,x - *),T(r,x - *) 



dr 



< Cp,K,Tt 



ap 



for some constant c P) k,t, where a > is the parameter in hypothesis (H5). 

Proof. We proceed as in the proof of |NQ07 Proposition 3.3]. For all n > 1, let ip n be 
the approximation of the identity defined in |NQ07 , Proof of Lemma 3.2]. Define the 
functions, G^ x (r,£) := (ip n (*) * {<Tij{u(t - r,*)) - a tj (u(t,x)))T(r,x - *))(£), T£(r,£) : = 
(ipn(*) * r(r, x — *))(£), and *^(r, £) = (ip n (*) * ^{f, x — *))(£)> where * is the convolution 
product and ^ is the measure defined in hypothesis (H5). Then, Gn X (r,*), T%(r, *), and 
tf£(r,*) belong to y(R d ;R). 
Then, for all n > 1, 



E 



(G^(r,*),r*(r,*))^dr 



JR d 



JR d 



G t 7 i x (r,y)f(y-z)T x n (r,z)dydzdr 



T(r,x- dy')(aij(u(t - r,y')) - aij(u(t,x)))T(r,x - dz') 



Using Minkowski's inequality with respect to the finite measure 





p- 


j dr 





4>n(y - y')f(y - z)ip n (z - z')dydz I dr, 



T(r, x — dy')T(r, x — dz') 
together with the Lipschitz property of a, we get that the last term is bounded by 

c !Kk I I I r^x-dy'^ElWuit-^y^-u^xW^r^x-dz') 

\ v 

V>n(y - y')f(y - z)^ n {z - z)dydz dr . 
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Therefore, by Lemma 13.21 we find that 



E 



< C P,K,T 



JR d 
X 



+ 



Cp,K,T 



r(r, x — dy')\\x — y'|| 74//2 r(r, x — dz') 
4>n{y ~ y')f(y - z)ip n (z - z')dydz)dr 
r(r, x — dy')r^^ 2 T{r, x — dz') 
Tpn(y - y')f(y - z)ip n (z - z')dydz]dr 
+ 



(^ n (r,*),T x n (r,*))^dr 



r^ 2 / \&T*(r)(d\ 2 »(dt;)dr 



Hence, by hypothesis (H5), we have that 



supE 

n>l 



< c p,K,T 



(G^(r,*),TUr,*))^dr 



($(r,*),r(r,*))jg>dr 



+ 



r 7/2 / \3?r(r)(0\ 2 n((%)dr 



< c P ,K,T(t aiP + e° 2P ) < c p , K)T e 



ap 



(3.4) 



Finally, since ^Vn(0 — ► 1 pointwise, appealing to Fatou's lemma, and (|3.4p . we conclude 
that 



E 



((Tij(u(t - r, *)) - <7ij(u(t,x)))r(r,x - *),T(r,x -*)) dr 



.-// 



< liminf E 

ri —>oo 

< c p , K , T e ap . 



(G^(r,*),r£(r, 



□ 



4 Existence and smoothness of the density 

In this section we prove Theorem 11.1( a) . 

First of all, we recall some elements of Malliavin calculus. Consider the Gaussian fam- 
ily {W(h),h S J%t} defined in Section 2, that is, a centered Gaussian process such that 
E[VF(/i)iy(g)] = (h,g),tf> T . Then, we can use the differential Malliavin calculus based on it 
(see for instance [N06]). We denote the Malliavin derivative by D = (D^ l \ D^'), which 
is an operator in L 2 (Q;J%r). For any m > 1, the domain of the iterated derivative D m in 
L p {n;je® m ) is denoted by B m >P, for any p > 2. We set D°° = n p >i n m >i B m ' p . Recall 
that for any differentiable random variable F and any r = (fx, ...,r m ) € [0, T] m , D m F{r) 
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is an element of ffl® m which will be denoted by D™F. We define the Malliavin matrix of 
F G (ID) 00 )™ by if = ((DFi, DFj)j& r )x<ij< m . 

In order to prove Theorem ll.l( a). we will use the following localized variant of Malliavin's 
absolute continuity theorem. 

Theorem 4.1. [BP98, Theorem 3.1] Let S m C M. k , m €E N, m > 1, be a sequence of open 
sets such that E m C S m +i and let F £ (D°°) fc such that for any q > 1 and m£N, 

E[(det lF )-n {F ^ m} ] <+oo. (4.1) 

Then the law of F admits a < ^°° density on the set £ = U m E m . 

The next result is the £>dimensional extension of (NQ07, Proposition 6.1]. Its proof 
follows exactly along the same lines working coordinate by coordinate, and is therefore 
ommitted. 

Proposition 4.2. Assume that hypotheses (HI) and (H2) are satisfied. Then, for every 
(t,x) € (0, T] x M. d , the random variable Ui(t,x) belongs to the space D°° ; for all i = 1, k. 
Moreover, the derivative Dui(t,x) = (D^Ui(t,x), D^Ui(t, x)) is an Jfffi -valued process 
that satisfies the following linear stochastic differential equation, for all j = 1, k: 

D^Ui(t,x) = (Tij(u(r, *))T(t - r,x - *) 



+ 



J^r(t - s,x - y)f^D^(a te (u(s,y))W\ds,dy) (42) 

F(t — s, dy)D^\bi(u(s, x — y))ds, 

for all r E [0, t], and is otherwise. Moreover, for all p > 1, m > 1 and i = 1, k, it holds 
that 



sup E 

(t,x)e{0,T]xR d 



\D m Ul (t,x)\\l^ k)i 



< +oo (4.3) 



Remark 4.3. Recall that the iterated derivative of Ui(t,x) satisfies also the fc-dimensional 
extension of equation |NQ07, (6.27)]. 



Recall that the stochastic integral on the right hand-side of (|4.2p must be understood by 
means of a Hilbert- valued stochastic integral (see |NQ07[ Section 3]) and the Hilbert-valued 
pathwise integral of (|4.2p is defined in |NQ07[ Section 5]. 

We next prove Theorem II .ll f a) and Remark II .21 
Proof of Theorem a). For each m £ N, m > 1, define the open set 

S m = L G R k : ||e T • a(y)\\ 2 > — , V£ G M fc , = l}. (4.4) 

By Proposition 14.21 it suffices to prove that condition (14. ip of Theorem 14.11 holds true in 
each S m . Indeed, oberve that S = U m S m = {deter ^ 0}. 



12 



Let (t,x) G]0,T] x M rf and m > 1 be fixed. Let c := |c, where c is the constant 
in hypothesis (H6). It suffices to prove that there exists 5o(m,c) > such that for all 
< 5 < 5q , and all p > 1 , 



(4.5) 



for some A > 0, and for some constant C > not depending on 5. This implies (14,1ft . taking 
p = f + 1 in (|43D . 



We write 



det7 u(ti!r ) > nnf S6R &.||^|| =1 (C r 7«(t )a; )0l • 



Let £ G M fc with |j£|| = 1, and fix e G (0, 1]. The inequality 

2^2 2 

11° + b\\jg>k > ^IMI.j^fc ~~ 2|[6||^fc, 



(4.6) 



(4.7) 



together with (|4.2p . gives 



> 



t 

t-e 



^2,D r (v,i(t,x))£i 



i=l 



dr > -M 
~ 3 



where 



i=l ^i" 6 

i=i 



(r,t,x,*) = ^2 / r(t-s,x-y)D r (o- tf (u(s,y))^(ds,dy) 




D r (bi(u(s, x - y))T(t - s, y)dyds. 

Now, assume that u(t,x) G S m . Similarly as in the proof of Lemma 13.31 define the 
smooth functions T n (r, £) := (i/j n (*) *T(r, *))(£)■ Then, adding and subtracting the term 
(£ T • a(u(t, x)))f into afi, and using inequality (|4.7p . we get that ^ > — 2=2^2, where 



4,1 = lim ^(^■a(«(t,xM 



1=1 



"1,2 



|^T n (r)(0|V(de)rfr, 

JO JR d 

lim / / / / T(t — r, x — dy')T(t — r,x — dz') 

x ( (C T • a(n(r,y'))) 4 ^ T ■ " (f ■ ) * 
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Note that we have added and subtracted a "local" term to make the ellipticity property 
appear (see (|4.4|) ). A similar idea is used in [MS99J for the stochastic wave equation in 
dimension 2. 

Then, using the fact that u(t,x) G T, m and hypothesis (H6), we get that 



1 



1 



M,i >~ I \J?T(r)(0\ 2 n(d0dr > c-e\ 
m J J R d m 

Now we find out upper bounds for the p-th moment, p > 1, of the terms and i^. 
We start treating =2^2- We write 

(£ T • <r(u(r,y)))i (£ T ■ <r(u(r,z)))i - (£ T ■ a(u(t,x))) 2 
= (C T ■ a(u(r, y)))i • (a(u(r, z)) - <r(u(t, x))j) . 

+ (£ T • a(u(t, x)))i (£ T ■ (a(u(r, y)) - a(u(t, x))))i. 

Then, proceeding as in the proof of Lemma 13,31 using the Lipschitz property of the coeffi- 
cients of a, and (|3,1|) . it yields that 



E 



SUp 1^*1,2 1 
?6R d :||e||=l 



We next treat s&i- Using the Cauchy-Schwarz inequality, for any p > 1, it yields that 

k r p-i 



E 



SUp \J$ 2 \ 



<c p ^T(E 

1=1 ^ 



<c p ^E / ||aj(r,t,a:,*)||;^ fc eir 

i=l L Jt-e 



\Vi(r,t,x, *)WL, k dr 



+ E 



/ \\Wi(r,t,x,*)\\ 2 ^ k dr 
Jo 



where 



V±(r,t,x,*) :=vf / r 

„ , Jt-r JR d 



t 



(t- s,x- y)D t - r (aa(u(s, y))W e (ds, dy), 



Wi(r,t,x,*) := / r(t - s,y)D t - r {bi(u(s,x - y))dyds. 

Jt~r JR d 

Now, using Holder's inequality, the boundedness of the coefficients of the derivatives of a, 
|NQ07t (3.13), (5.26)], and hypothesis CLU}, we get that 



E 



\Vi(r,t,x, *)\\%> h dr 



<c v ^- x sup E[\\D^. Ui (t-s,y)\\% ck ] 



(s,y)6(0,e)xl 



xi// |^r( s )(OlV(de)rf« 

10 JR d 



< c p e 2 ^ 2P . 
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Moreover, using Holder's inequality, the boundedness of the partial derivatives of the coef- 
ficients of b, hypotheses (|1.4|) and (|1.6|) . and NQj07j (5.17), (5.26)], for the second term in 
(|4.8I) corresponding to the Hilbert- valued pathwise integral, we obtain that 

pi 



E 



\Wi(r,t,x,*)\\jg, k dr 



<c p e p 1 sup E[||A— «t(* - s,'i 

(s,i/)e(0,e)xR d 



x|l / T(s,y)dyds 
10 JR d 



< c.e (72+1)p . 



Hence, we conclude that, for any p > 1, 



E 



sup 



< c p e ap . 



Appealing to H4.6I) . we have proved that, on the set {u(t, x) G S m }, 

k 



det7u(t,a) > 



rn 



(4.9) 



where I is a random variable such that for all p > 1, E[|I| P ] < e ap , and 0,77 > are the 
parameter of hypotheses (H5) and (H6). 

We now choose e = e(5, m, c) in such a way that 5 l / k = ^e' 7 , an d since e < 1, 5 < 5q := 
(2m)-fc, xhen using (|4.9p . we conclude that for all < 5 < 5 and p > 1, 

= C{m,c,p)5 Xp , 



Pi det7 u(ti:r )l {u ( ti:c)eEm } < 5 \ < 



JL e V _ gl/k 



with A = (?7 < a). This proves (|4.5p . 



□ 



Proof of Remark \1. 6 A If a is strongly elliptic, we can use the same computations above 
without localizing on the set T, m , and then apply [DKN09a, Proposition 3.5] with Z := 
inf m=1 (£ T "/ u(ttX )0, *l,e = Y 2 , e = sup||£|| =1 (|^ 2 | + Ki,2|), eo = 1, at = a 2 = V, and 
Pi = fh = ol (rj < a), to conclude that for any p > 1, 



E[(det 7u 



(t,aO, 



where the constant c Pt T,K is clearly uniform over (i, x) in any nondegenerate compact rectan- 
gle A x 5 C]0, T] x M d . This, together with Proposition EO] and |DKN09al Theorem 3.1 and 
Proposition 3.4], proves that under the strong ellipticity condition, for any (t, x) G]0, T] xM. d , 
the law of the random vector u(t, x) admits a < ^°° density which is uniformly bounded over 
any nondegenerate compact rectangle Ax B c]0,T] x M. d . □ 



5 Criterium for the strict positivity of the density 

Given T > 0, a predictable processes g G J^r, h 6 Jf k , and z G R fc , we define the process 
W = (W 1 , W k ) as 

Wi(l m h j ) = W j (l m h j )+z j [ (hjW^is^jeds, j = l,...,k,te[0,T}. 

Jo 
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We set W t (h) = £* =1 W j (l [0it] hj). Then, by [DZ921 Theorem 10.14], {W t ,t G [0,T]} is a 
cylindrical Wiener process in J$? k on the probability space (f2, J^,P), where 



— (w) = j z , to & n, 



where 

J z = exp 



J>/ / gj (s,y)W^ds,dy)-lj2z] [ \\c?(s, *)\\%, ds). (5.1) 
« io it" 2 'H Jo J 



Then, for any predictable process Z G L 2 (f2 x [0,T];J^' k ) and j = 1, fc, it yields that 

f T I Z j (s,y)W j (ds,dy) = C f Z j (s,y)W^ds,dy) + Zj C \z j (s,*),g^s,*))^ds. 
JO Jwi d JO JR d Jo 

Now let u z (t,x) the solution to equation (jl.2p with respect to the cylindrical Wiener 
process W, that is, for i = l, k, 

k 



uf(t,x)=J2[ [ T{t-s,x-y)a ij {u z {s,y))W ] {ds,dy) 
~{ Jo JR d 

k ft 

+ 5> / < r (* 

3=1 J ° 



3=1 
k 

s,x--)a l:j (u z ( S ,*)),g3(s,*))^ ( 5 - 2 ) 



+ / / bi(u z (t - s,x - y))T(s,dy)ds. 



Then, the law of u under P coincides with the law of u z under P, that is, for all i = 1, k, 

E[ui(t,x)] = E[uf(t,x)J z ]. (5.3) 

Given a sequence {g n }n>i of predictable processes in J#r, h G J^ k , and z G R k , let 
Un(t, x) be the solution to equation (jl.2p with respect to the cylindrical Wiener process 
{WT,t€ [0,T]}, where W?(h) = £ -=i ^(l^fy), and 

^ nj (l[o,t]^) = ^'(l[o,t]^) +*J [\h,(*),gi(s,*))^ds. 

Jo 

Set the k x k matrix (p^(t,x) := d z ii^(t, x), and the Hessian matrix of the random vector 
u*(t, a;), ip^t^x) := d1u^{t,x) (which is a tensor of order 3). We denote by || • || the norm 
of a n x n matrix A defined as 

\\A\\ = sup ||A£||. 

,11*11=1 



We next provide the main result of this section which is a criterium for the strict 
positivity of a density, which was proved in [BP981 Theorem 3.3] for the case where = 
L 2 ([0, 1]; M fc ). Our case follows similarly along the same lines as theirs, but for the sake of 
completeness we provide its proof. This criterium uses the following quantitative version of 
the classical inverse function theorem. Let B(x;r) denote the ball of M k with center x and 
radius r > 0. 
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Lemma 5.1. [BP98, Lemma 3.2] and [N98J Lemma 4.2.1] For any /3, 5 > 0, there exists 
constants R G]0, 1], a > 0, such that any ff 2 mapping g : M. k i— > R fc satisfying 

|det</(0)|>- and sup (|<?(x)| + \g'(x)\ + \g"(x)\) < fi, 

P \\x\\<8 

is a diffeomorphism from a neighborhood of contained in the ball B(0;R) onto the ball 
B(g(0);a). 

Theorem 5.2. Assume that for all (t,x) G]0,T] x IR^, that the law of random vector u(t, x) 
has a continuous density p t:X on X an open subset o/R fe , that is, 

lz(y)P °u~ 1 (t,x)(dy) = lY,{y)pt,x{y)dy. 

Assume that for all (t, x) G]0, T] xM. d , there exists a sequence of predictable processes {<7 n }n>i 
in J%t, and positive constants c\, C2, r$ and 5 such that, for all y G Supp(P o u _1 (i, x)) H S ; 

(i) limsup^^pi (\\u(t,x) -y\\ < r) n (det<p° n (t,x) > ci) I > 0, for all r G]0,r ]. 



(ii) lim^ooP |sup|| z ||<,5(||^(t,x)|| + \\^(t, x)\\) < c 2 | =1. 
Then pt tX is a strictly positive function in X. 

Proof. Fix (t,x) G]0,T] x R d and y G Supp(P o u^fax)) n S, and set 
A n := i(||«(t,x) -y|| < r) n (dety>°(t,z) > cy) 



n ^ sup^(||^(t,a;)|| + |K(t,x)||) <c 2 J j. 

Then hypotheses (i) and (ii) imply that there exists ro > such that for each r g]0, ro], 
there exists n G N such that 

P{An} > 0. (5.4) 
Let r and n be fixed such that ()5.4p holds and note that, on A n , 

sup ||u*(i,a;)|| < + Sc 2 < \\y\\ + r + 5c 2 := k. 

\\*\\<8 

Then by Lemma [5.1|, there exists a > depending on c\,C2,S,k and R g]0, 1], such that for 
all u G A n the mapping z i— > u^(t,x) is a diffeomorphism between an open neighborhood 
V^(u;) of in M fc contained in the ball B(0;R), and the ball B(u(t,x);a). We choose r 
sufficiently small such that r < a, and i? small enough such that u G A n and z G V^(u;) 
imply that det^(i, x) > y. 
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Let g : R fc i— > R be a strictly positive continuous function such that JJg*. g(z)dz = 1 and 
let / : R fc nKbea nonnegative continuous and bounded function with support contained 
in E. Then by (|5.3p . we have that 

E[f(u(t,x)}= [ g(z)E[f(u z n (t,x))J z ]dz 



!a„ / g(z)f(u z n (t,x))J z dz 



> E 



> E 



U„ / g(z)f{u z n {t,x))J z dz 



Then, using the theorem of change of variables, we get that 



E[/(n(t,x)] >E 



1a 



> 



B{u{t,x);a) |det<£; 

f(v)O n (v)dv. 



——rg((u v n (t,x)) 1 )J {<(t ^ )) -idv 

n\ l i x )\ 



where 



(v) :=E 



and </3 : R + i— > [0, 1] is a continuous function such that 1[ 0)7 .] < 93 < l[o,a]- 

Note that from (|5.4p and the fact that x) —y\\ < r} C A n , we have that n (y) > 0. 

Moreover, because the function 



v i-» tp{\\u(t,x) -w||)inff i det w f 1 )" J (^(*,a:))-i ) 1 

is a.s. continuous and bounded by 1, from Lebesgue's dominated convergence theorem, we 
conclude that n is continuous. 

Hence, for each y 6 Supp(Poit -1 (t, sc))nS, we have proved the existence of a continuous 
function On : R fc 1— ► R + such that #n(y) > 0, and for all / : R fc 1— > R nonnegative continuous 
and bounded function with support contained in S, 



f(z)p t , x (z)dz > / f{z)0l(z)dz. 



(5.6) 



This implies that pt, x (y) > 0. On the other hand, as n (y) > and the function n is 
continuous, there exists a neighborhood U y of y where n does not vanish. Hence, we 
deduce from (|5.6p that for all y G Supp(P o u (t, x)) n £ 

c Supp(P o u^it.x)) n s, 

which implies that Supp(P o a;)) is an open set, and thus contains E. This implies 

that for all y £ £, pt, x (y) > 0, and the theorem is proved. □ 
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6 Proof of Theorem 11.1( b) 



The aim of this section is to prove Theorem 11.1( b). For this, we assume that T is a 
nonnegative real valued function such that hypotheses (jl.4fl . (|1.3p . and (Hi), i=2,...,6, 
hold. We will prove that conditions (i) and (ii) of Theorem 15.21 are satisfied. 

Proof of Theorem M.lV b). Let (t, x) £]0, T] x W 1 be fixed. Consider the sequence of pre- 
dictable processes {g n }n>i in J&T, defined by 



We are going to prove that assumptions (i) and (ii) of Theorem 15.21 are satisfied for this 
sequence {g n } and for any y € Supp(P o u~ l (t,x)) n X, where £ = {detu ^ 0}. This will 
prove the theorem. 

Proof of (i). Using the interpretation of the Malliavin derivative as a directional derivative 
we have that 



J o 

Then by (|5.2p and the stochastic differential equation satisfied by the derivative (|4.2p . we 
obtain that 



a; 



i(s,z) = c n 1 l [t _ 2 -n j4 ](s)r(i - s,x- z), n > 1, j = 1, k, 



(6.1) 



where 





where 



r* tifj (t,x) 




Jo 




\ F(t- s,x - y) d k o-u(u z n (s,y)) 

JO \Jr JR d £k=1 





X ( A 0) «,fc0, * - V)),9n(r, *))jeT(t - s, dy)ds dr. 



Note that gh,{r,*) = ^[t-2~ n ,t]( r )9n(r,*), and that D T {u z n j(s,y)) = if s < r. Hence, using 
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this facts and Fubini's theorem, it yields that 

k 



^n,i,j(t, x ) = I I T(t-s,x-y) Y] d k a ie (u z n (s, y)) 

Jt-2-» JR d i k=l 

x (D® (u z n>k (s, y)),gl(r, *))^dr ) W n > e (ds, dy), 

ft p _k 
/n,i,j{^ x ) = Y] d k bi{u z n {s, x - y)) 

Jt~2-" JR d k=1 

x ( I {D$Xu\ * k (s,x - y)),gi(r,*))^dr)r(t - s,dy)ds. 



Therefore, we have proved that 

l Pn,i,jV'i x ) ^n.i.j \Pt x ) ~r" ^n,i,jV'i x ) ^n,i,j\^i x ) c^n,i,jy J i x )i 

where 



^n,i,j(t, x ) ■= c n 1 / (ff -(t^(t-r 1 *))r(r,x-*),r(r,x-*)> jr dr, 
J o 

ft f k 

^n,i,j(ti x ) ■= / / T(t-s,x-y) V d k a u {u z n (s, y))<p z k j (s,y)W i (ds,dy), 

k ,-t 

^n,i,j{^ x ) '■= Y] HI ( r (t - s,x - *)d k a u (u^(s,*))<p z nk j(s,*),g e n (s,*))^ds, 
e >k =i Jt z 

ft f 

</n,i,j(M) = / / y2d k bi(u z n (s,x-y))^ z nAj (s,x-y)T(t- s,dy)ds. 

Jt-2- n JU d k=1 

Next we study upper bounds for the p- moments, p > 1, of the last three terms. Ap- 
pealing to |NQ07j (3.11)], (|Q|> . and hypothesis (fTB]) . and using the fact that the partial 
derivatives of the coefficients of a are bounded, we get that for all p > 1, 

k 

Efl^MH < c p ]T sup E[|^ fcJ («,y)H 
k=1 (s,y)e[0,T]xRd 

x(f [ \^T{t-s,x-*)(0\ 2 KdOds) P 

\Jt-2- n JR d J 
k 

< Cp , T 2-^Yl su p nwi^^n 

k=1 (s,yM0,T]xR d 

We now use Minkowski's inequality with respect to the finite measure T(s, dy)ds, the 
boundedness of the partial derivatives of the coefficients of b, and hypothesis (II. 4jl . to see 
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that 



E[!A«(^)n<c P E SU P V[\<ti,kj(*,v)\ P ]( f [ r(s,dy)ds 

k=1 (s,y)e[0,T]xR d \J0 JR d 

k 

<c PiT 2- n Pj2 sup E[\rt Aj (s,y)n 

fe=1 ( S)2 ,)e[0,T]xR« 

In order to bound the third term jtf^i j(A x )-> assume that \\z\\ < 5 for some 5 > 0. Then, 
using the Cauchy-Schwarz inequality, and the fact that the derivatives of a are bounded, it 
yields that 



E[K M (t,x)n<c p , T ^^E 



k=l 



t-2~ n 



l ( ^n,fc,j(s ) *)r(t ~S,X- *)\\%ds 



Then, proceeding as in the proof of Lemma 13.31 and using hypothesis (II. 3D . we get that 



k=1 (s,y)e[0,T]xR d 



X 
k 



2 -n 

JR d 



\P/2 

\&T(s)(0\ 2 KdH)ds\ 



A . =1 (a,!/)e[0,r]xR<' 



Hence, we have proved that for all p > 1, 

k k 



^EM tkjj (t,x)\r] < c p ]TE[|^ fcj (t,z)n + c ; iT (2-™^/ 2 + 2~ np + F) 



k=l 



k=l 



xV sup E[|^. fcJ (s,y)| ? 

fc=1 (s,j/)e[0,T]xRd 



(6.2) 



On the other hand, using the Cauchy-Schwarz inequality, and proceeding as in the proof 
of Lemma l3.3l using the boundedness of the coefficients of a, and hypothesis (|1.3p . we obtain 
that 



E[|^-(t,x)H<E 



Wij«(t -r,*))T(r,x- *)\\ 2 ^dr 



(6.3) 



We next choose n sufficiently large and 5 sufficiently small such that for any p > 1 

c' PjT {2~ n ^ 2 + 2~ np + 5 P ) < i (6.4) 
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Then, (|C2j) . (|6T3D . and ([631) im p!y that 

V sup sup E[\(p z nkJ (s,y)\ p ] < Cp, T . 

fc=1 p||<(5(s,y)e[0,T]xR d 

Moreover, the calculations above also show that 



(6.5) 



E 



In particular, when z = 0, this implies that for all p > 1 



<c p , T (2-^ +2-"P + 5P). 



E 



Ki, J -(f,z)-<i, J -(M)l p 



< c PjT (2- n72p/2 + 2" 



(6.6) 



On the other hand, we write 

A) 



W n iAt,%) = aij(u(t,x)) + nii j{t,x), 



where 



2 -n 



{<Jij{u(t - r, *)) - CTjj(ii(t,x)))r(r,x - *),T(r,x - *) ) dr. 



'o \ / je 

By Lemma 13.31 and hypothesis (H6), we conclude that for all p > 1, 



E \0 n ,ij(t,x)\* <c PtK , T 2-^ a -^ (a> V ). (6.7) 

Now, since y £ Supp(P o ii _1 (i, x)) H X, there exists ro > such that for all < r < ro, 
B(y; r) C S, and P {u(t, x) G r)} > 0. 
Moreover, det a(y) > (say). Hence, for all < r < ro, 

P { (\\u(t, x) - y\\ < r) n (det a(u(t, x)) > 2d) I > 0, 



(6.8) 



where 



ci := - inf a(z). 

2 zeS(y;r) 



Furthermore, we have that 
det (p^ (t, x) > det a(u(t, x)) — det ^a(u(t, x)) + R n (t, x)\ — det a(u(t, x)) 
where R n {t,x) is the k x k matrix R n (t,x) = & n (t,x) + ip^(t,x) — W®(t,x), and 
det ya(u(t, x)) + R n (t, x)\ — det a(u(t, x)) 

= det a(u(t, x))< det I I& + <r~ (u(t, x))R n (t, x) ) — 1 



(6.9) 



(6.10) 
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where 1^ denotes the k x k identity matrix. 
Moreover, observe that (for k > 2), 



det flk + o 1 ('u(t, x))R n (t, x)\ — 1 = det <r a;)) det R n (t, x) 

fc-l fc 

+ J] ^(u^a!))^^!))^ •••(a- 1 (n(t,x)) J R n (t,x)) 

fc=l h,...,i k =l,ii<---<i k 



(6.11) 



Therefore, by the properties of the determinants (|6.10p and ([6. lip , and using (|6.6p and 
(|6.7p , we conclude that for all p > 1 , 



lim E 

n— >oo 



0. 



det ia(u(t, x)) + R n (t, x) J — det a(u(t, x)) 
Finally, (|H^]l . {Eg) , and (l6~T2l) imply that 

limsupP i (||n(t,x) - y|| < r) n (det</>°(£,x) > c x ) \ > 0, 

which proves (i). 

Proof of {li). We start proving that there exist c > and 5 > 0, such that 

lim P J sup \\v z n (t,x)\\ <c\ = 1. 
n ^°° l|MI<5 J 

Let i>^(t,:c) denote the matrix such that each entry is solution to the affine equation 



(6.12) 



(6.13) 



£,k=l ^° 

(6.14) 

Define 

Ii,j(f,w) = Y] zi / (r(t - s,z - *)5 fc <T rf «(s,*))/ fcj (s,*),^(s,*))^ds. 

Then, for all it; £ $7 and ||z|| < 5, using the boundedness of the partial derivatives of the 
coefficients a, the definition of the g n 's, and using the approximation argument of the proof 
of Lemma 13.31 an d hypothesis (|1.3p . it is easily shown that, for all j = 1, k, 



IW'^lloo ^ Sc T\\f: 



3 Noo- 



Hence, for 5 small enough, fj t— > Ij(f,u) is a contraction for all uj£!1. Consequently, for 5 
small enough, equation (|6.14|) has a unique solution such that, by the boundedness of the 
coefficients of a and its derivatives, for all n > 1 and l|z|| < 5, it holds that 



\W(t,x)\\<K'\\W*(t,x)\\<K. 
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Therefore, for all p > 1, 



P<j sup \\<ti(t,x)\\>2K <P sup |K(i,x)-<(i,x)|| >*T 

||z||<<5 J U\z\\<S 



(6.15) 



Now, using the equation for the random field (p^(t,x) — v^(t,x) and similar computations 
as in the proof of (i), it yields that, for all p > 1 and ||z|| < 5, 

E[\\tp* n (t,x)-v z n (t,xW] <c p , T (2-^/ 2 + 2-^) sup E[M(s,yW) 

(s,y)e[0,T]xR d 

+ c' p 5 p sup E[\\<p*(s,y)-v*(s,y)n 

(s,y)e[0,T]xK d 

Choose 5 small enough such that c'5 p < |, and use ()6.5[) to see that 

sup sup E[||<( S ,y) - <( S)2 /)|H < c PiT ,k(2-^ 2P / 2 + 2-"f). 

IMI<<5 O,j/)e[0,T]xR d 



(6.16) 



Observe that this implies that 



lim sup P <^ \\<p z n {t,x)\\ > 2K \ = 0. 

n ^°°\\z\\<5 



Therefore, in order to prove (|6.13|) one needs to check the uniformity in z. Let z and 
z' such that ||z|| V \\z'\\ < 5. Then, using the same calculations that lead to (|6.2p . the 
Lipschitz property of a' and b' , and subtracting and adding the terms S~*(t,x), J^(t,x), 
and ^ z (t,x) with u z n (s,y) and (p z n {t,x), we get that 

E [\\ri(t,x) - f^(t,x)|H < c p E [||^(t,x) - ^(t,x)H + c; T)if (2-»™/ 2 + 2"»f + P) 



x( sup E[||<( S ,y)-<( S ,y)|| p ] 

,(s,j/)e[0,T]xR d 



+ sup E 

(s,y)e[0,T]xR d 



\u z n (s,y) -u z n (s,y))\\ p \\ip z n (s, 



Choosing 5 and n as in (|6.4|) . applying the Cauchy-Schwarz inequality, and appealing to 
(I6.5|) . it yields that 

sup E y) - (s, y) f] < c p sup E [||^( S , y) - *tf (s, y) f] 

(s,y)e[0,T]xR d (s,y)e[0,T]xR d 

, ,V2 

+ Cp,T,K SUp E 

(s,j/)e[0,T]xR d 

Proceeding along the same lines, we get that 



\<{s,y) -u z n(s,y)\\ 2p 



sup E[\\v z n (s,y)-v z n (s,y)\\»] < c p sup E[\\W n z (s,y) - W z (s,yW] 

(s,?/)e[0,T]xR d (s,y)e[0,T]xR d 

1/2 



+ cj>,T,K sup E 

(s,y)e[0,T]xR d 



-«n(s,y)ll 2p 



24 



We now claim that for all p > 1, 

sup E[\\e n (s,y)-u^(s,y)\\P}<c p , T , K \\z-z'\\ p - (6-17) 

(s,y)e[0,T]xR d 

Indeed, using the Lipschitz property of the coefficients of a and b, together with [NQ07, 
(3.9), (5.15)], we obtain that 

E[\\u* n (t,x)-tf(t,xW]<c p \\z-z'\\P 

+ c p ,t,k [ su P E[|K( S ,y)-<'( S ,y)|H f \&T(t - s)(0\ 2 ^)ds 
Jo j/eM d ^R d 

+ c p ,t,k sup E [\\u*(t - s,y) - u„(t - s,y)\\ p ] / T(s,dy)ds. 

JO y£R d JR d 

Thus, hypotheses (|1.4|) and (|1.3p . and Gronwall's lemma prove (|6.17|) . 

We next use the approximation argument of Lemma 13.31 together with the Lipschitz 
property of a, and (|6.17|) . to show that 

sup E[\\W*(s,y)-W*'(s,yW}<c PiT7 K\\z-z , \\i>. 

(s,i/)e[0,T]xM d 

Therefore, we have proved that 

E[||<(t,x)-<'(t,x)|H+E[||<(t,x)-<'(t,x)|H<c p , Tii ,||z-z'|r, 
which, together with (16.161) . shows that for all p > 1, 



lim E 

n^oo 



sup \\(f%(t,x) -v*(t,x)\\ p 

\\z\\<5 



0. 



This, together with (16. 151) conclude the proof of (I6.13|) . 

The proof of (ii) for ip^(t,x) follows along the same lines, therefore we only give the 
main steps. Let 

d 2 

Then, using the interpretation of the Malliavin matrix as a directional derivative, we have 
that ^ ^ 

Vn,i,j, m {^x) = / / {D^D^(u z nyi (t,x)),gi(r,*) ® g™(s,*)).^^drds, 
Jo Jo 

where the Jif® ^-valued process (u^ j(t, x)), satisfies the following linear stochas- 

tic differential equation 

D^D^(u z n:i (t,x)) =F(t-r,x- *)0< m W«( r >*))) + T(t-s,x- *)£>«> (a im (u(s, *))) 
+ [ [ r(t-w,x-y)^Di m ~>D^(a ie (u(w,y)))W e (dw,dy) 

JrVs JM. d £=1 

+ 11 T{t-w,dy)Di m ^D^{b i {u(w,x-y)))dw. 



rVs 
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Using the chain rule and the stochastic differential equation satisfied by the first derivative, 
one can compute the different terms of ip z n ; • m (t, x) as we did for ip* i .(t, x). Also, a similar 
but simpler equation holds for d z v*(t,x). 
One first shows that for all p > 1, 



lim sup E 

n ^°°|! 2 ||<«5 



l/ n (t,x)-d z v z n (t,xW 



0. 



Then one estimates the pth-moments of the differences ip^(t, x) — ^ (t, x) and d z v^(t, x) — 
d z >v* {t-> x ) as we did for (p^(t,x) in order to get the desired result. □ 



7 Proof of Theorem 11.3 



The existence and smoothness of the density follows from |NQ07 Theorem 6.2]. Hence, we 
only need to prove the strict positivity. For this, one applies Theorem 15.21 as in Theorem 
11.1( b) taking S = R. In this case, in order to prove hypothesis (i) of Theorem 15.21 one 
proceeds as in Theorem 11.1( b). using hypotheses (HI), (H2), and (II. 6p . to show that for 
all p > 1, 



lim E 

n— >oo 



0. 



We next find a lower bound for W®(t,x). For all m > 1, define the familly of smooth 
functions T^(r,e) := ftM*)*r(r, *))(£) and 4f(r,0 := (^m(*)*«r(«(t-r, *))T(r, x-*))(0 
as in the proof of Lemma 13.31 Then, using the non-degeneracy assumption on a, we get 
that 



W n °(t,x) 



lim c" 



1 / <^*(r,*),r*(r, *)>.**■ 



o 



-i 



lim c 

m^oo / q 



> c. 

Hence, this implies that 

<p n (t,x)>c-\<p n (t,x)-#*{t,x)\. 
Finally, these assertions imply that for all y G Supp(P o n _1 (t,x)) 

limsupP \(\u(t,x)-y\ < r) n (y>° (t, x) > £)) > 0, 

which proves (i). The proof of (ii) follows exactly as in Theorem ll.lf b). 



8 Examples 

In this Section we apply our Theorems 11.11 and 11.31 to spatially homogeneous non-linear 
stochastic heat and wave equations. 
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8.1 The stochastic heat equation 

Let r be the fundamental solution of the deterministic heat equation in W d with null initial 
conditions, that is, 

Ar = o. 

dt 2 

Then, T is given by the Gaussian density 

T(t,x) = (2iTt)- d/2 expf 



2t 
and 

^r(t)(e) = ex P (-2^neii 2 ). 

Note that hypothesis (II, 4p is clearly satisfied. At is pointed out in |NQ07 Example 4.3], 
hypothesis (|1.3|) holds if and only if 

/ 7TTWT <+ »' ^ 

jR d I 1 + U\\ ) 

Moreover, it was proved in [SS00J that under the stronger condition: there exists e G (0, 1) 
such that 

(TTIFF < + °°' <8 - 2) 

hypotheses (|1.5j) and (jl.6j) are satisfied for all 71 6 (0, 1 — e] and 72 G (0, 1— e], and hypothesis 
(H4) holds for all 74 G (0,2(1 - e)]. Condition ([T7]) is clearly satisfied for all 73 G (0,1]. 
Also, condition (H6) is satisfied for any r/ > 1 (see |NQ07, Remark 6.3]). 

Observe that under condition (|8,2p . by Lemma 13.21 the trajectories of the solution u of 
equation (|1.1|) where L is the heat operator are 7/2-Holder continuous in time, and 7-Holder 
continuous in space, for all 7 G (0, 1 — e]. 

We finally study hypothesis (H5). For this, we assume that the spatially homogeneous 
covariance is defined as the Riesz kernel f(x) = for < (3 < (2 A d). That is, 

fj-(d^) = \\£,\\P~ d d!;. In this case, the integrability condition (|8.2|) holds for all e > 

For all h G [0, T], using the changing variables [y = ^yjf, z = 3 ^=r], it yields that 

fh 

\\y — z\\~^Y(r, y)T(r, z) dydzdr 

d 

r -fi/2 f f ||~ _ z\\-Pr(l,y)T(l,z)dydzdr = ch^ . 

jR d jR d 

2-/3 
2 • 

Moreover, by [SS00] . hypothesis (|1.5j) holds for all 71 G (0, ^5^), and hypothesis (H4) 
holds for all 74 G (0, 2 — (3). Condition (jl.7p is clearly satisfied for 73 = 1. Hence, 7 = 71. 

Furthermore, hypothesis (H5) holds with a% = + ^ and 02 = + -j-- Indeed, 
using the change of variables [y = z = -^], for all r G [0, 1], we have that 

z \\^ \\u ~ z ll _ ^r(r, z)T(r, y)dzdydr 

f f „~,i2±u~ — ~x^,-. -x , 2-|3 1 74 



JR d 



Thus, hypotheses (jl.6p and (H6) are satisfied for 72 = rj 



JR d 



r 1 I I \\--\\ II//— i|| ^r(l, z)T(l, y)dzdydr = ct -' ! 
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and 



/ r 2 1 / / ||y — z\\ ^r(r, z)Y(r,y)dzdydr 

JO JR d JR d 

= / r ^ _ 2 / / ||y - z\\~ l3 Y(l,z)Y(l,y)dzdydr = ct^ +1 2 . 
JO JR d JR d 

Observe that a x A a 2 = ^ + < |(2 - /3) < (2 - (5) = (2 72 ) A ( 72 + 1). 

Hence, applying Theorem 1 we have the following result. 

Theorem 8.1. Let u be the solution of equation (jl.ip . where L is the heat operator and 
the spatial covariance of the noise is given by a Riesz kernel. Assume that the coefficients 
of the equation satisfy hypothesis (H2). Then, for all (t,x) €]0,T] x M. d , the law of random 
vector u(t,x) admits a strictly positive density on the set Y := {deter ^ 0}. 

Finally, as a consequence of Theorem ll.3l we prove the strict positivity of the density of 
a single stochastic heat equation in any space dimension, under the integrability condition 
(USD- 
Theorem 8.2. Let u be the solution of equation (jl.ip with k = 1, and L the heat operator. 
Assume that a and b are c &°° bounded functions with bounded derivatives, and \o~\ > c > 0. 
Assume also that (18.20 holds. Then, for all (t,x) g]0, T] x IR^, the law of random variable 
u(t, x) admits a c £°° strictly positive density on M. 



8.2 The stochastic wave equation 



Let Yd be the fundamental solution of the deterministic wave equation in R with null initial 
conditions, that is, 



d 2 T n 



dt 2 

Then, for d = 1,2,3, Yd is given, respectively, by 

T\(t,x) = -l{|x|<t}- 

T 2 (t,x) = ^-\\x\\% 1/2 . 

T 3 (t,x) = ±- t a t , 

where at denotes the surface measure on the three-dimensional sphere of radius t. Note 
that only for d € {1,2,3}, Yd defines a nonnegative measure. Moreover, for all d > 1, 

*r,w(a - sm(2 ' t|lfll) 



2tKII 

Note that hypothesis (jl.4p is clearly satisfied. At is pointed out in |NQ07[ Example 4.2], 
hypothesis (jl.3p holds if and only if (|8.ip is satisfied. 
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Moreover, if we assume condition (|8.2p . then hypothesis (jl.5p holds for all 71 G (0, 2(1 — 
e)] (see [ggOOl Lemma 3]), hypothesis ([IS]) holds for all 72 G (0,3 - 2e] (see |qS04| Lemma 



3]), and condition (H4) is satisfied for all 74 G (0, 2(1 — e)] (see [SS00, Lemma 3]). Finally, 
condition CC?]) holds for 73 = 2 (see |QS04j (A.4)]) (H6) is satisfied for rj = 3 (see (QSM 
(A.3)]). 

In particular under condition (|8.2p . by Lemma 13.21 the trajectories of the solution u of 
equation 111. IB where L is the wave operator with d G {1,2,3} are 7-Holder continuous in 
time and in space, for all 7 G (0, 1 — e]. 

As in the heat operator case, in order to check hypothesis (H5), we assume that the 
spatially homogeneous covariance is defined as the Riesz kernel. 

Observe that in this case, as 71 G (0, 2 — f3), 72 G (0, 3-/3), and 73 = 2, we have that 
7 G (0,2-/3). Moreover, 74 G (0,2-/3). 

Notice also that 

\^(r)(i)\ 2 Uf~ d d£dr<T^ T I \^T d (r)(0\ 2 U\f~ d d£dr < +00. 



>0 JM. d JO JR d 

Moreover, changing variables, for all r G [0, 1] and d G {1, 2, 3}, we have that 
(*(r,*),r d (r,*)).x.dr= f #T d {r)(Z) ^ffjP ll^f ~" <%dr 



Jo 

*), r d (l, *))jr [ T r 2 "^ dr = cr 3 -^ 
Jo 



and 







ri / \^T d (r)(t)\ 2 Uf- d d£dr = cT 3 -P + 2. 



Hence, hypothesis (H5) holds taking a = 3 — (3 + 2 ^ i (< 4 — 1/3), and choosing 72 such 
that 4 — |/3 < (272) A (72 + 1). Notice that this is possible because 72 G (0, 3 — f3). 

Finally, in order to have a = 3 — /3 + > 3 = rj (required in hypothesis (H6)), it is 
necessary that 7 A 74 > 2/3. But as 7 A 74 < 2 — (3, this is only possible if /3 G (0, |). 

Hence, appealing to Theorem ll.il we have proved the following result. 

Theorem 8.3. Let u be the solution of equation (II. ip . where L is the wave operator and the 
spatial covariance of the noise is given by a Riesz kernel with (3 G (0, |) and d G {1,2,3}. 
Assume that the coefficients of the equation satisfy hypothesis (H2). Then, for all (t,x) G 
]0, T] x M. d , the law of random vector u(t,x) admits a strictly positive density on the 
set T := {deter / 0}. 



Finally, as a consequence of Theorem 11.31 we have the following result for the case of a 
single wave equation. 

Theorem 8.4. Let u be the solution of equation (II. ip with k = 1, and L the wave operator 
with spatial dimension d G {1,2,3}. Assume that a and b are c to co bounded functions 
with bounded derivatives, and \o~\ > c > 0. Assume also that (|8.2|) holds. Then, for all 
(t,x) g]0, T] x M. d , the law of random variable u(t,x) admits a < ^°° strictly positive density 
on R. 
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